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Abstract 

A  1/32°  global  ocean  nowcast/forecast  system  has  been  developed  by  the  Naval  Research  Laboratory  at  the  Stennis  Space  Center. 
It  started  running  at  the  Naval  Oceanographic  Office  in  near  real-time  on  1  Nov.  2003  and  has  been  running  daily  in  real-time 
since  1  Mar  2005,  It  became  an  operational  system  on  6  March  2006,  replacing  the  existing  1/16°  system  which  ceased  operation  on 
1 2  March  2006.  Both  systems  use  the  NRL  Layered  Ocean  Model  (NLOM)  with  assimilation  of  sea  surface  height  from  satellite 
altimeters  and  sea  surface  temperature  from  multi-channel  satellite  infrared  radiometers.  Real-time  and  archived  results  are  available 
online  at  http ;//www, ocean. nrlsscmavy.mil/globaL_nlom.  The  1/32°  system  has  improvements  over  the  earlier  system  that  can  be 
grouped  into  two  categories;  (1)  better  resolution  and  representation  of  dynamical  processes  and  (2)  design  modifications.  The  design 
modifications  arc  the  result  of  accrued  knowledge  since  the  development  of  the  earlier  1/16°  system.  The  improved  horizontal 
resolution  of  the  1/32°  system  has  significant  dynamical  benefits  which  increase  the  ability  of  the  model  to  accurately  nowcast  and 
skillfully  forecast.  At  the  finer  resolution,  current  pathways  and  their  transports  become  more  accurate,  the  sea  surface  height  (SSH) 
variability  increases  and  becomes  more  realistic  and  even  the  global  ocean  circulation  experiences  some  changes  (including  inter¬ 
basin  exehange).  These  improvements  make  the  1  /32°  system  a  better  dynamical  interpolator  of  assimilated  satellite  altimeter  track 
data,  using  a  one-day  model  forecast  as  the  first  guess.  The  result  is  quantitatively  more  accurate  nowcasts,  as  is  illustrated  by  several 
model-data  comparisons.  Based  on  comparisons  with  ocean  color  imagery  in  the  northwestern  Arabian  Sea  and  the  Gulf  of  Oman,  the 
1/32°  system  has  even  demonstrated  the  ability  to  map  small  eddies,  25-75  km  in  diameter,  with  70%  reliability  and  a  median  eddy 
center  location  error  of  22,5  km,  a  surprising  and  unanticipated  result  from  assimilation  of  altimeter  track  data.  For  all  of  the  eddies 
(50%  small  eddies),  the  reliability  was  80%  and  the  median  eddy  center  location  error  was  29  km.  The  1/32°  system  also  exhibits 
improved  forecast  skill  in  relation  to  the  1/16°  system.  This  is  due  to  (a)  a  more  accurate  initial  condition  for  the  forecast  and  (ft)  better 
resolution  and  representation  of  critical  dynamical  processes  (such  as  upper  ocean  -  topographic  coupling  via  mesoscale  flow 
instabilities)  which  allow  the  model  to  more  accurately  evolve  these  features  in  time  while  running  in  forecast  mode  (forecast 
atmospheric  forcing  for  the  first  5  days,  then  gradually  reverting  toward  climatology  for  the  remainder  of  the  30-day  forecast  period). 
At  1/32°  resolution,  forecast  SSH  generally  compares  better  with  unassimilated  observations  and  the  anomaly  correlation  of  the 
forecast  SSH  exceeds  that  from  persistence  by  a  larger  amount  than  found  in  the  1/16°  system. 
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1*  Introduction 

The  development  of  a  daia-asstmilative  1/32°  global 
ocean  nowcast/forecast  system  is  the  culmination  of  a 
long-standing  effort,  which  includes  the  associated  basic 
research  and  exploratory  development.  For  operational 
purposes,  the  U.S.  Navy  niche  in  ocean  modeling  and 
prediction  is  high  horizontal  resolution.  For  global  and 
basin-scales  that  means,  for  example,  the  ocean  models 
must  have  sufficient  resolution  to  (1)  depict  narrow 
oceanic  current  systems  and  their  meandering,  (2)  rep* 
resent  sharp  ocean  fronts  and  eddies,  (3)  represent  geo¬ 
metric  and  topographic  features  such  as  narrow  passages 
and  islands,  (4)  simulate  upper  ocean  -  topographic  cou¬ 
pling  via  flow  instabilities  that  is  strong  enough  to  ac¬ 
curately  depict  the  pathways  of  many  current  systems  in 
the  world  ocean  and  (5)  provide  boundary  conditions  for 
even  higher  resolution  coastal  models. 

How  much  resolution  is  enough?  Dramatic  improve¬ 
ments  are  found  as  the  grid  resolution  is  increased  up  to 
1/32°  (^  3.5  km  at  mid- latitudes),  but  only  modest 
changes  are  found  when  the  resolution  is  increased  to 
1/64°  (Hurlburt  and  Hogan,  2000).  A  similar  result  was 
found  in  the  Japan/East  Sea  (Hogan  and  Hurlburt,  2000), 
As  the  resolution  is  increased,  the  pathways  of  current 
systems  become  much  more  accurate,  ocean  fronts  be¬ 
come  sharper  and  extend  much  farther  to  the  east  across 
major  ocean  basins,  eddies  fill  the  entire  model  domain, 
the  large-scale  shape  of  basin-scale  ocean  gyres  is  al¬ 
tered,  and  the  global  ocean  circulation  is  changed,  in¬ 
cluding  inter-basin  exchanges,  and  the  pathways  of 
global  scale  currents  and  their  transports.  One  of  the  key 
aspects  of  model  dynamics  that  drives  the  resolution 
requirement  is  the  need  to  resolve  the  vorticity  dynamics 
of  baroctinic  instability,  not  just  the  resulting  mesoscale 
eddies.  This  is  particularly  evident  when  studying  the 
impact  of  ocean  model  resolution  on  upper  ocean  - 
topographic  coupling  via  mesoscale  flow  instabilities. 
Baroclinie  instability  is  very  effective  in  transferring 
energy  from  the  upper  ocean  to  the  abyssal  ocean  and  is  a 
critical  mechanism  in  driving  many  deep  mean  flows  (as 
well  as  time  varying  flows)  (Holland  and  Lin,  1975; 
Holland,  1978;  Rhines  and  Holland,  1979;  and  other 
references  in  this  paragraph).  In  ocean  models  the  energy 
transfer  and  resulting  abyssal  currents  increase  greatly  up 
to  -  1/32°  resolution.  In  turn  abyssal  currents  at  mid 
and  high  latitudes  can  be  very  effective  in  steering  upper 
ocean  currents  and  improving  simulations  of  their  path¬ 
ways  (Hurlburt  and  Thompson,  1980;  Hurlburt  et  al., 
1996;  Hurlburt  and  Metzger,  1998;  Hogan  and  Hurlburt, 
2000;  Tilburg  et  al.,  2001;  Hogan  and  Hurlburt,  2005), 
From  all  of  the  preceding  discussion,  the  1/32°  model 


should  have  improved  dynamical  interpolation  skill 
when  assimilating  data,  and  it  should  demonstrate  im¬ 
proved  forecast  skill.  However,  the  preceding  results 
indicate  that  increases  in  grid  resolution  beyond  1/32° 
will  yield  only  modest  improvements  in  most  deep-water 
regions  of  the  World  Ocean  (much  higher  resolution  is 
needed  in  coastal  regions).  Hence,  our  target  resolution 
globally  is  1/32°.  In  relation  to  the  size  of  ocean  eddies, 
3/32°  (~  3,5  km  mid-latitude)  resolution  for  the  ocean  is 
similar  to  the  resolution  currently  used  by  most  weather 
forecasting  models  in  relation  to  the  size  of  atmospheric 
highs  and  lows,  and  in  relation  to  the  corresponding  first 
internal  mode  radii  of  deformation  for  the  ocean  and  the 
atmosphere. 

Applications  for  ocean  nowcast/forecast  systems  in¬ 
clude  assimilation  and  synthesis  of  global  satellite  sur¬ 
face  data;  ocean  prediction;  optimum  track  ship  routing; 
search  and  rescue;  anti-submarine  warfare  and  surveil¬ 
lance;  tactical  planning;  high  resolution  boundary  con¬ 
ditions  that  are  essential  for  even  higher  resolution 
coastal  models;  sea  surface  temperature  for  long  range 
weather  prediction;  inputs  to  ice  models,  biochemical 
models  and  shipboard  environmental  products;  environ¬ 
mental  simulation  and  synthetic  environments;  observ¬ 
ing  system  simulation  and  assessment;  ocean  research 
and  education;  pollution  and  tracer  tracking;  fisheries 
and  marine  resource  management;  impact  on  ocean 
structures  such  as  oil  rigs,  and  inputs  to  water  quality 
assessment,  Data-assimil alive  eddy-resolving  ocean 
prediction  systems  can  be  very  effective  in  using  real¬ 
time  satellite  altimeter  data  to  map  and  forecast  the 
“ocean  weather"  (e.g.  ocean  eddies  and  the  meandering 
currents  and  fronts)  as  demonstrated  in  earlier  studies 
(e.g,  Hurlburt  et  al„  2000;  Smedstad  et  al„  2003).  Sea 
surface  temperature  is  assimilated  as  well. 

Here  we  discuss  the  operational  1/32°  near-global 
NLOM  nowcast/forecast  system  and  the  value  added  of 
the  resolution  increase  over  the  earlier  l/l 6a  system 
discussed  in  Smedstad  et  al,  (2003),  The  1/16^  system 
has  been  running  continuously  in  real-time  at  the  Naval 
Oceanographic  Office  (NAVOCEANO),  Stennis  Space 
Center,  MS,  since  18  Oct.  2000  and  it  became  an 
operational  system  on  27  Sept,  2001.  The  1/32°  system 
has  been  running  in  near  real-time  since  1  Nov.  2003  and 
daily  in  real-time  since  l  Mar.  2005.  It  became  an  op¬ 
erational  system  on  6  March  2006,  replacing  the  existing 
1/16°  system  which  ceased  operation  on  12  March  2006, 
Results  from  both  systems  are  available  online  at  http:// 
www.ocean.nrlssc.navy.mil/globaLnlom.  Both  systems 
use  the  NRL  Layered  Ocean  Model  (NLOM)  with 
assimilation  of  sea  surface  height  from  satellite  alti¬ 
meters  and  sea  surface  temperature  from  multi-channel 
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satellite  infrared  radiometers.  The  next  section  gives  an 
overview  of  this  system.  Improvements  in  the  nowcast 
accuracy  and  forecast  skill  of  the  1/32°  over  the  1/16° 
system  are  discussed  in  Section  3,  and  the  conclusions 
are  summarized  in  Section  4, 

2.  The  nowcast/ forecast  system 

2.  L  The  NRL  layered  ocean  model  (NLOM) 

The  model  component  of  the  ocean  prediction  system 
is  based  on  the  primitive  equation  model  of  Hurlburt  and 
Thompson  (1980)  but  with  greatly  expanded  capability 
(Wallcraft,  1991;  Wallcraft  and  Moore,  1997;  Moore  and 
Wallcraft,  1998;  Wallcraft  et  al,  2003).  The  model  has  a 
nearly  global  domain  that  extends  from  72°S  to  65°N. 
The  horizontal  resolution  of  each  model  variable  is  1/32° 
in  latitude  by  45/1024°  in  longitude  or  —  3.5  km  at  mid¬ 
latitudes,  which  is  eddy -resolving.  It  has  6  dynamical 
layers  plus  the  mixed  layer  and  vertically  compressed  but 
otherwise  realistic  bottom  topography  that  is  confined  to 
the  lowest  layer  of  the  model,  except  that  flows  through 
straits  are  constrained  to  small  values  below  the  sill 
depth.  The  model  has  lateral  boundaries  that  follow  the 
200-m  isobath  with  a  few  exceptions,  such  as  the  shallow 
straits  around  the  Japan/East  Sea.  At  the  solid  boundaries 
kinematic  and  no-slip  boundary  conditions  are  used. 
Much  of  the  deep-water  formation  in  the  far  North 
Atlantic  is  parameterized  via  observational ly  based  flows 
through  northern  boundary  ports  at  the  Davis  Strait  and 
the  three  straits  between  southern  Greenland  and 
Scotland. 

In  the  model,  the  prognostic  variables  are  layer  den¬ 
sity,  layer  thickness,  layer  volume  transport  per  unit 
width  (layer  velocity  times  layer  thickness),  sea  surface 
temperature  (SST)  and  mixed  layer  depth  (MLD).  The 
model  has  a  free  surface,  corresponding  to  sea  surface 
height  (SSH)  (i,e.T  a  variable  observed  by  satellite  alti¬ 
metry).  The  Kraus-Tumer  type  mixed  layer  is  not 
confined  within  the  upper  dynamical  layer.  To  some 
extent  it  is  independent  of  the  dynamical  layers,  but  is  not 
entirely  passive.  In  particular,  (1)  a  deep  mixed  layer  can 
distribute  surface  forcing  across  multiple  dynamical 
layers,  (2)  thermal  expansion  is  based  on  the  mixed  layer 
temperature  Tm  rather  than  layer  1  temperature,  and  (3) 
surface  heat  flux  depends  on  All  three  factors  can 
change  the  steric  sea  surface  height  anomaly.  The 
embedded  mixed  layer  model  used  here  carries  prog¬ 
nostic  equations  for  SST  and  MLD  and  is  discussed  in 
detail  by  Wallcraft  et  ah  (2003). 

Below  the  mixed  layer,  the  density  of  the  top  5  dy¬ 
namical  layers  is  relaxed  toward  the  annual  mean 


climatological  density  of  that  layer  except  for  layer  1 , 
which  is  relaxed  toward  a  monthly  mean  climatology 
interpolated  to  daily  values  because  of  the  significant 
seasonal  cycle  within  that  layer.  The  Modular  Ocean 
Data  Assimilation  System  (MODAS)  climatology  (Fox 
et  al.,  2002)  is  used  for  this  purpose.  Unlike  other  models 
with  fixed  levels  in  the  vertical,  such  relaxation  does  not 
significantly  damp  the  anomalies  because  in  NLOM 
most  of  the  information  about  circulation  anomalies  is 
carried  by  layer  thickness  variations,  not  density  vari¬ 
ations.  For  example,  NLOM  maintained  a  Rossby  wave 
generated  by  the  1982-83  El  Nino  for  at  least  a  decade 
(Jacobs  et  al,  1994)  without  oceanic  data  assimilation 
except  for  the  relaxation  to  climatological  density  within 
layers. 

Additional  details  concerning  NLOM  can  be  found  in 
Hurlburt  et  al.  ( 1 996),  Metzger  and  Hurlburt  (1996),  and 
Shriver  and  Hurlburt  (1997).  Shriver  and  Hurlburt 
(1997)  discuss  the  ability  of  global  NLOM  to  simulate 
the  global  overturning  circulation,  including  the  cross- 
interfacial  mixing  scheme  based  on  oxygen  saturation. 
Mixing  also  occurs  when  there  is  model  layer  outcrop¬ 
ping.  Kara  and  Hurlburt  (2006)  discuss  NLOM  SST 
simulation  skill  (with  no  assimilation  of  SST  data)  on 
time-scales  ranging  from  daily  mean  values  to  inter- 
annual 

2.2,  The  assimilation  scheme 

With  a  model  grid  the  size  of  the  1/32°  NLOM 
(8192x4608  per  layer),  it  is  essential  to  use  an  assim¬ 
ilation  technique  that  will  not  result  in  a  large  increase  in 
the  model  run  time.  The  scheme  used  is  similar  to  the 
incremental  updating  technique  described  in  Smedstad 
and  Fox  ( 1 994)  and  is  the  same  data  assimilation  scheme 
used  in  the  current  operational  1/16°  NLOM  system 
(Smedstad  et  al,  2003).  This  technique  increases  the 
model  run  time  by  about  50%  compared  to  the  model 
running  without  assimilation.  More  sophisticated  tech¬ 
niques  like  the  ensemble  Kalman  filter  (Even sen,  1994) 
may  increase  the  run  time  of  the  assimilation  system  by  a 
factor  of  up  to  1 00* 

The  SSH  data  that  are  assimilated  into  the  model  are 
satellite  altimeter  SSH  deviations  from  a  temporal  mean 
plus  the  mean  SSH  described  in  the  next  subsection. 
Currently  Jason- 1,  Geo$at-Follow-On  (GFO)  and  Envi- 
sat  are  used  in  the  real-time  system.  However,  data  from 
ERS-2,  GFO  and  Topex/Poseidon  were  used  for  the  June 
2001  -May  2002  hindcast  results  discussed  in  this  paper 
In  the  data  assimilation  scheme  the  model  SSH  is  used 
as  a  first  guess  in  an  OI  deviation  analysis  with  a  3-day 
data  window.  Anisotropic,  spatially  varying  mesoscale 
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covariance  functions  determined  from  altimeter  data 
(Jacobs  et  ah,  200 1 )  are  used  in  the  01  deviation  analysis. 
The  surface  information  is  transferred  to  the  deeper 
layers  via  the  statistical  inference  technique  of  Hurlburt 
et  al.  ( 1 99G).  This  technique  updates  the  pressure  field  in 
all  layers  below  the  surface.  When  the  coefficients  tor  the 
statistical  inference  were  calculated,  care  was  taken  to 
use  a  recent  time  window  of  model  fields,  namely  1 997- 
2002  at  a  15-day  interval.  A  geostrophic  velocity  cor¬ 
rection  is  calculated  from  the  pressure  changes.  The 
velocity  correction  is  not  performed  between  5°  latitude 
and  the  equator.  Between  5  and  8°  the  correction  is 
gradually  increased  to  full  strength  using  a  hyperbolic 
tangent  function.  To  further  reduce  the  creation  of  grav¬ 
ity  waves,  these  corrections  are  then  used  incrementally 
to  update  the  model  variables. 

The  SST  assimilation  consists  of  a  relaxation  of  model 
SST  (mixed  layer  temperature,  Tm)  toward  the  MODAS 
MCSST  analysis  (T0).  MODAS  is  an  operational  system 
running  daily  at  NAVOCEANO  (Fox  et  al.,  2002).  The 
relaxation  process  is  treated  as  a  heat  flux  where  the  flux 
depends  on  an  e-folding  time-scale  (r^-3  h  during  the 
assimilation)  and  the  mixed  layer  depth  (/rm).  In  addition, 
the  total  heat  flux  (Qy  from  relaxation  and  the  atmosphere 
(0a))  is  limited  to  0mm  and  Qmax  of  ±  1 000  WrrT2,  Le. 

Or  =  max  (0mm.  minima,,  0A  +  Cpphmr^{T0-Tm))), 

where  Cp  is  the  specific  heat  of  water  and  p  is  water 
density.  A  benefit  of  treating  the  SST  assimilation  as  a  heat 
flux  (with  a  flux  limiter)  is  demonstrated  in  Zamudio  et  al, 
(2002).  In  that  paper  1/16°  global  NLOM  nowcasts  show 
cold  SST  beneath  Hurricane  Juliette  despite  warm 
MODAS  SSTs  derived  from  the  surrounding  region  due 
to  lack  of  data  beneath  the  cloud  cover  of  the  hurricane. 

During  each  30-day  forecast,  the  forecast  SST  is 
relaxed  toward  climatologically-corrected  persistence  of 
the  nowcast.  In  this  case  the  relaxation  time-scale  is  1/4 
the  elapsed  forecast  length  (e.g.,  1  week  for  a  4-week 
forecast).  This  procedure  is  used  because  the  atmospher¬ 
ic  forecast  is  only  available  out  to  5  days.  During  the 
forecast  the  model  makes  the  SST  look  more  fluidic  and 
helps  keep  SST  fronts  better  aligned  with  currents  and 
SSH  fronts  than  climatologically-corrected  persistence. 
It  also  helps  keep  eddies  in  SSH  and  SST  in  phase. 

2.3.  The  mean  sea  surface  height 

The  mean  sea  surface  height  used  in  the  assimilation 
scheme  (Fig.  t)  was  developed  using  the  same  pro¬ 
cedures  as  in  Smedstad  et  at.  (2003).  It  is  important  to 
have  a  mean  which  has  the  major  ocean  currents  and 


associated  SSH  fronts  more  sharply  defined  than  is  fea¬ 
sible  from  hydrographic  climatologies.  As  a  starting 
point,  the  mean  SSH  from  the  atmospherically-forced 
model  (with  no  ocean  data  assimilation  except  for  the 
climatological  density  relaxation)  was  calculated  over 
the  time  period  of  the  statistical  inference  statistics 
( 1 997-2002).  The  accuracy  of  this  mean  is  substantially 
enhanced  by  the  climatological  density  relaxation.  The 
mean  was  then  compared  to  available  sources  of  in¬ 
formation,  including  mean  surface  dynamic  height  over 
the  time  period  of  the  satellite  observations,  mean  frontal 
locations  determined  from  satellite  IR,  mean  SSH  de¬ 
termined  from  nearly  simultaneous  AXBT  underflights 
of  altimeter  tracks  (Mitchell  et  aL,  1990)  and  SSH 
variability  from  satellite  altimeter  data.  The  model  mean 
was  then  modified  via  a  rubber  sheeting  technique  and 
other  methods  designed  to  operate  on  SSH  fields  (Carnes 
et  ah,  1996).  The  Gulf  Stream,  Gulf  of  Mexico,  Kuroshio, 
and  Japan/East  Sea  were  the  only  regions  where  mod¬ 
ifications  were  made. 

2.4.  The  wind  and  thermal  forcing 

The  wind  and  thermal  forcing  used  in  the  1/32°  model 
are  the  same  as  used  in  the  previous  1/16°  operational 
system  (Smedstad  et  al.,  2003).  The  wind  stress  used  is  a 
hybrid  of  the  Fleet  Numerical  Meteorology  and  Oceano¬ 
graphy  Center  (FNMOC)  Navy  Operational  Global  Atmo¬ 
spheric  Prediction  System  (NOGAPS)  (Hogan  and 
Rosmond,  1991;  RosmondetaL,  2002)  and  the  Hellerman 
and  Rosenstein  (HR)  (1983).  The  long  term  mean 
(August  1990- July  1999  for  FNMOC)  is  subtracted 
from  the  FNMOC  stresses  and  replaced  by  the  annual 
mean  from  H  R,  From  the  atmospheric  model  we  get  the 
forcing  up  to  a  5-day  forecast.  The  forcing  for  the  30-day 
forecast  gradually  reverts  to  climatology  beyond  5  days. 
The  last  forecast  record  is  weighted  with  the  contempo¬ 
raneous  climatological  values  with  the  weight  of  the 
forecast  decreasing  to  zero  over  a  1 0-day  time  span.  In 
Section  3,5,3  we  discuss  the  effect  of  this  on  ocean 
forecasts.  The  thermal  forcing  also  comes  from  FNMOC 
NOGAPS.  The  latent  and  sensible  heat  flux  formulation 
is  replaced  by  that  of  Kara  et  ah  (2002),  which  makes  use 
of  model  SST. 

Kara  et  al.  (2003)  and  Kara  and  Hurlburt  (2006) 
describe  experiments  showing  the  ability  of  NLOM  to 
accurately  predict  SST  over  most  of  the  global  ocean 
without  relaxation  to  or  assimilation  of  SST  data.  These 
experiments  show  that  NLOM  predicts  annual  mean  SST 
with  an  RMS  error  of  <  0.5  °C  for  the  nearly  global  model 
domain  in  comparison  with  the  Comprehensive  Ocean- 
Atmosphere  Data  Set  (CO ADS)  SST  climatology.  Using 
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Fig.  I.  Slightly  modified  mean  SSH  used  in  the  assimilation  of  satellite  altimeter  SSH  (discussed  in  Section  2.2)  in  the  (A)  1/32°  and  (B)  1/16“ 
systems.  The  1/32*  mean  was  computed  over  the  time  period  1997-2002  and  the  1/16°  mean  was  computed  over  the  period  1993-1999.  This  SSH 
field  is  added  to  the  altimeter  SSH  deviations  from  a  temporal  mean. 

inter-annual  simulations  with  6-hourly  forcing,  Kara  and 
Hurlburt  (2006)  found  that  a  1/8°  version  of  the  global 
model  gave  a  median  RMS  SST  error  of  .6  °C  in  com¬ 
parison  with  194-year  long  daily  time  series  (1990- 
1998)  measured  at  moored  buoys  in  the  Pacific  from  the 
National  Data  Buoy  Center  and  the  Tropical  Atmosphere 
Ocean  (TAO)  array  (McPhaden,  1995), 

2,5.  The  altimeter  and  SST  data 

The  real-time  altimeter  data  assimilated  into  the 
model  are  delivered  via  N  AVOCEANO’s  Altimeter  Data 
Fusion  Center  Currently,  data  from  the  GFO,  Jason- 1 
and  Envisat  altimeters  are  available  within  48  h.  Better 
orbit  corrections  are  available  with  a  slightly  longer 
delay.  To  pick  up  data  with  improved  orbits  and  other 
delayed  data,  the  real-time  assimilation  cycle  restarts 
5  days  prior  to  the  nowcast  time.  However,  many  of  the 
model -data  comparisons  and  the  comparisons  between 
the  1/32°  and  1/16°  global  ocean  prediction  systems  use 


the  period  June  200 1 -May  2002  when  altimeter  data 
from  ERS-2,  GFO  and  Topex/Poseidon  were  available. 
For  SST,  the  model  assimilates  MOD  AS  SST  analyses  of 
satellite  MCSST  data  (in  situ  SST  data  are  not  used). 
These  OI  analyses  are  a  daily  operational  product  of 
NAVOCEANO  (Fox  et  ah,  2002). 

2.6,  The  pre-operational  cycle 

The  operational  nowcast/forecast  system  is  currently 
running  daily  on  316  IBM  SP4+  processors  at  the 
NAVOCEANO  Major  Shared  Resource  Center  (MSRC), 
Each  week,  the  Wednesday  nowcast  is  used  to  initialize  a 
30-day  model  forecast  with  the  atmospheric  forcing  de¬ 
scribed  in  Section  2,4. 

3.  Results  and  discussion 

The  improvements  noted  in  the  1/32°  compared  to  the 
1/16°  system  are  attributable  to  better  resolution  and 


a 
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representation  of  dynamical  processes,  with  some  im¬ 
provement  due  to  design  modifications  as  well.  The 
design  modifications  are  the  result  of  accrued  knowledge 
since  the  development  of  the  earlier  1/16°  system.  To 
illustrate  the  value  added  of  the  1/32°  system  over  the 
earlier  1/16°  system,  several  model-data  and  model- 
model  comparisons  are  discussed  in  this  section. 

3.1.  Arabian  sea  SSH  and  currents  (with  and  without 
data  assimilation)  ur  Sea  WiFS  ocean  color 

Satellite  altimetry  provides  the  key  available  data 
type  that  allows  a  data-assimiiative,  eddy-resolving 
ocean  model  to  accurately  map  the  mesoscale  (e.g., 
Hurlburt,  1984;  Hurlburt  et  ah,  2000).  Here,  this  is 
illustrated  in  the  northern  Arabian  Sea  and  the  Gulf  of 
Oman  by  comparing  eddies  seen  in  three  ocean  model 
nowcasts  with  eddies  seen  in  a  SeaWIFS  ocean  color 
image  (Fig.  2).  Such  a  cloud-free  image  during  a  very 
large  outbreak  of  high  chlorophyll  is  an  unusual  event  in 
this  region.  In  the  model  with  atmospheric  forcing  but 
without  assimilation  of  the  SSH  and  SST  data,  the 
position  of  the  eddies  is  poorly  represented.  This  is 
expected  because  these  features  are  largely  non-deter- 
ministic  due  to  flow  instabilities.  With  the  assimilation, 
the  positions  of  eddies  at  both  I/!  6°  and  1/32° 
resolution  are  in  good  agreement  with  the  advective 
effects  in  the  SeaWiFS  imagery,  e.g.  note  the  column  of 
4  alternating  counterclockwise  and  clockwise  eddies 
extending  southward  from  the  coast  of  Iran.  These 
eddies  were  depicted  by  assimilating  GFO,  ERS-2  and 
Jason- 1  altimeter  data  into  the  1/32°  ocean  model  but 
only  data  from  GFO  and  ERS-2  in  the  operational  1/16° 
model,  since  Jason- 1  data  were  not  yet  available  in  real¬ 
time  at  NAVOCEANO.  By  using  the  finer  resolution 
model  with  data  from  a  third  altimeter,  some  of  these 
features  and  some  smaller  features  are  more  accurately 
positioned  and  shaped,  e.g.,  the  eddy  centered  near 
1 8°N,  62°E  and  an  observed  cyclonic  eddy  in  the  Gulf 
of  Oman  (at  approximately  57.5°E,  25°N),  which  is 
anti-cyclonic  in  the  1/16°  model. 

3.1.1.  A  quantitative  analysis  of  model  eddy  center 
location  error 

A  quantitative  analysis  of  model  eddy  center  location 
error  in  comparison  to  the  ocean  color  is  provided  in 
Fig.  3A  and  Table  L  All  but  one  eddy  with  a  clearly 
defined  center  in  the  ocean  color  were  used.  The  ex* 
ception  is  the  cyclonic  swirl  located  inside  a  large  anti* 
cyclonic  eddy  (northeast  part  of  #3).  The  first  search  for 
eddies  in  the  ocean  color  map  was  done  independently  of 
the  ocean  model  results.  This  search  identified  19  eddies 


in  the  ocean  color.  As  a  second  pass,  the  1/32°  and  1/16° 
data-assimiiative  ocean  models  were  used  as  an  aid  in 
finding  more  eddies  in  the  ocean  color.  This  search 
yielded  one  more  eddy,  the  smallest  one  (#20).  Eddy  #20 
is  about  1/4°  in  diameter  in  both  the  ocean  color  and  the 
data-assimiiative  1732°  model  (the  only  one  to  capture 
it).  Some  eddies,  especially  anti-cyclonic  eddies,  do  not 
have  a  clearly  defined  center  in  the  ocean  color.  Only 
35%  of  the  eddies  identified  in  the  ocean  color  were  anti- 
cyclonic.  Anti-cyclonic  eddies  generally  lack  the  up- 
welling  and  inward  spiral  commonly  associated  with 
cyclonic  eddies,  although  rapidly  strengthening  anti- 
cyclonic  eddies  can  have  anti-cyclonic  inflow  (Hurlburt 
and  Thompson,  1976)  and  sinking.  Sinking  is  a  char¬ 
acteristic  unfavorable  to  phytoplankton. 

The  latitudes  and  longitudes  of  the  eddy  centers  in  the 
ocean  color  were  determined  first  and  without  any  ref¬ 
erence  to  the  model  results.  The  eddy  centers  in  the  model 
results  were  determined  primarily  from  the  velocity  field, 
secondarily  using  the  SSH  field  because  the  ocean  color 
field  is  strongly  influenced  by  advection.  The  model  eddy 
centers  were  determined  as  independently  as  possible 
from  the  ocean  color  and  each  other.  In  this  process  the 
model  and  ocean  color  eddies  were  compared  only  to  the 
extent  needed  to  match  them  up  geographically  and 
determine  that  they  had  the  same  sense  of  rotation  and 
similar  diameters.  To  qualify  as  a  match  the  model  eddies 
had  to  meet  all  three  of  these  criteria.  To  match 
geographically  the  model  eddy  and  the  ocean  color 
eddy  had  to  overlap  or  be  capable  of  overlapping  on  a 
week-long  time-scale  using  a  differential  propagation 
speed  up  to  10  km/day.  The  preceding  is  based  on  the 
assumption  that  some  eddies  were  depicted  by  assimila¬ 
tion  of  altimeter  data  at  least  a  wreek  in  the  past.  For  the 
smaller  eddies,  this  gives  an  upper  bound  on  the  eddy 
center  location  error  of-  100  km.  However,  55  km  is  the 
largest  error  in  Table  I  for  any  small  eddy  where  this 
upper  bound  value  applies.  The  accuracy  in  determining 
the  locations  of  the  eddy  centers  is  typically  10-15  km  for 
both  the  ocean  color  and  the  models. 

The  eddies  in  Table  1  are  listed  in  the  order  of  de¬ 
creasing  diameter  as  seen  in  the  ocean  color.  The  purpose 
is  to  aid  in  assessing  the  eddy  mapping  skill  of  the  ocean 
models  assimilating  satellite  altimeter  SSH  as  a  function 
of  decreasing  eddy  size.  The  first  four  eddies  are  -  2°  or 
more  in  diameter,  eddies  5-8—1°,  12-16-1/2°  and 
1 7-20—  1/4°.  With  a  few  exceptions,  model  eddies  1 2- 
19  are  typically  -  1/2°  to  3/4°  in  diameter  with  the 
remainder  generally  about  the  same  as  depicted  in  the 
ocean  color.  A  1/4°  model  eddy  like  #18  or  #20  is  8  grid 
points  in  latitude  by  6  points  in  longitude  in  the  1/32° 
model,  about  the  smallest  it  could  realistically  depict. 
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Fig.  2.  A  comparison  of  (A)  chlorophyll  concentration  from  SeaWiFS  (2-6  Oct  2002  latest  cloud  free  pixel  composite  with  most  data  from  6  Oct) 
with  (B-D)  NLOM  SSH  (with  surface  currents  overlain)  from  6  Oct  2002,  (B)  1/32°  model  with  SSH  and  SST  assimilation,  (C)  l/l  6°  model  w  ith  the 
assimilation,  and  (D)  t/32c  model  without  the  assimilation.  All  were  subject  to  atmospheric  forcing.  In  (D)  the  positions  of  most  eddies  are  poorly 
represented  because  they  are  predominantly  non^deterministic  due  to  flow  instabilities.  With  the  assimilation  of  SSH  data  from  satellite  altimeters,  the 
eddy  positions  show  good  agreement  with  adveedve  effects  in  SeaWiFS  imagery,  e.g,  note  die  column  of  4  alternating  counterclockwise  and 
clockwise  eddies  extending  southward  from  the  Iran  coast.  In  order  to  facilitate  inter-comparison,  a  constant  offset  (7  cm)  was  subtracted  from  (C)  to 
better  reference  it  to  the  areal  mean  in  (B)  and  (D). 


The  1/32°  NLOM  with  assimilation  of  SSH  from  3 
altimeters  was  able  to  depict  90%  of  the  20  eddies  seen 
in  the  ocean  color,  all  of  the  largest  10,  8  of  the  10 
smallest  vs  70%  depicted  by  the  1/16°  NLOM  with 
assimilation  from  2  altimeters,  8  of  the  10  largest  and  6 
of  the  10  smallest.  The  1/32°  NLOM  with  no  synoptic 
ocean  data  assimilation  depicted  35%  of  the  eddies,  2  of 
the  10  largest  and  5  of  the  10  smallest.  The  median 
center  location  error  for  the  eddies  depicted  is  29  km  for 


the  1/32°  NLOM  with  assimilation,  35.5  km  for  the  1/ 
16°  NLOM  with  assimilation  and  42  km  for  the  1/32° 
NLOM  with  no  assimilation.  All  20  of  the  eddies  seen  in 
the  ocean  color  were  depicted  by  at  least  one  of  the  two 
data-assimilative  models.  The  1/32°  NLOM  with 
assimilation  had  the  lowest  center  location  error  for 
57,5%  of  the  eddies,  the  1/16°  NLOM  for  27,5%  of  the 
eddies  and  the  1/32°  NLOM  with  no  assimilation  for 
15%  of  the  eddies. 
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Fig,  3.  (A)  Same  as  Fig.  2  A  with  clearly  defined  eddy  centers  numbered  in  older  of  decreasing  eddy  size  (see  Table  l )>  Number  color  varies  only  for 
visual  clarity.  (B,  C)  Same  as  Fig.  2 B,  C  with  the  most  recent  week  of  observed  altimeter  tracks  overlaid  as  solid  lines  and  the  remaining  tracks  as 
dashed  lines,  red  for  ERS-2,  black  forGFO,  and  white  for  Jason- 1.  At  that  time  Jason- 1  data  were  not  available  in  real-time  at  NAVOCEANO  and 
thus  not  assimilated  by  the  operational  1/16°  NLOM  system,  but  the  tracks  arc  overlaid  on  (C)  for  reference.  (D)  Same  as  (B)  but  for  1  week  earlier. 


3.1.2.  Effectiveness  of  altimeter  data  assimilated  by  the 
models  in  depicting  small  eddies ,  25-75  km  in  diameter 
How  effective  is  the  assimilation  of  satellite  altimeter 
data  in  depicting  the  smaller  eddies?  The  1/32°  NLOMn 
with  no  assimilation  of  synoptic  ocean  data  contains 
numerous  small  eddies  and  one  might  expect  that  some  of 
these  would  match  up  with  real  eddies  seen  in  the  ocean 
color  by  chance,  because  of  geometric  constraints  (e.g.  for 


eddies  near  the  boundary)  and/or  because  of  the  atmo¬ 
spheric  forcing.  The  ten  smallest  eddies  shown  in  the  ocean 
color  (1 1-20)  will  be  used  to  investigate  this  issue.  Only 
two  are  adjacent  to  a  coastal  boundary,  #14  and  #20.  Eddy 
#14  is  depicted  by  all  three  models  in  Table  1  and  with 
similar  location  error,  while  #20  is  only  present  in  1/32° 
NLOM  with  the  altimeter  data  assimilation.  This  suggests 
that  a  boundajy  indentation  in  combination  with  wind 
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Table  I 


Eddy  center  location  errors  in  ocean  prediction  models  compared  to 
ocean  color  from  SeaWiFS  in  the  northwestern  Arabian  Sea  and  Gulf 
of  Oman 


Ocean  color  A 
eddy  ID#  or 

C 

Ocean  color 
eddy  center 
location 

1/32° 

NLOM 

1/16° 

NLOM 

1/32* 

NLOMn  no 
assim 

°N 

°E 

Eddy  center  position  error  in  km 

I 

C 

23.55 

60.2 

18 

35 

NP 

2 

A 

18,3 

62,0 

28 

103 

NP 

3 

A 

21.15 

60.25 

58 

44 

NP 

4 

C 

19,65 

60.4 

12 

43 

45 

5 

C 

16,8 

55,65 

42 

35 

31 

6 

c 

17.0 

58.9 

17 

42 

NP 

7 

c 

16.7 

57.9 

79 

53 

NP 

8 

A 

18.0 

57.6 

40 

NP 

NP 

9 

c 

25.1 

57,6 

39 

NP 

NP 

10 

A 

24,7 

62.5 

35 

30 

NP 

n 

c 

23.7 

62.3 

22 

* 

42 

12 

c 

19.3 

58.8 

30 

30 

NP 

13 

A 

19.225 

59.35 

11 

35 

NP 

14 

A 

25,3 

58.55 

33 

36 

30 

15 

c 

24,1 
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26 
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12 

NP 

NP 

%  of  eddies  present 

90% 

70% 

35% 

Median  eddy  center  position  error  in  km 

29 

35.5 

42 

%  of  eddies  with  most  accurate  position 

57.5% 

27.5% 

15% 

Notes: 

The  ocean  color  eddy  ID  numbers  are  plotted  on  Fig.  2a,  Eddies  are 
listed  in  order  of  decreasing  size  as  depicted  by  the  ocean  color.  Eddy 
position  measurement  error  is  10- 15  km  in  both  the  ocean  color  and 
the  models. 

A,  C:  A  for  anti -cyclonic  eddies*  C  for  cyclonic  eddies, 

I  /32°  NLOM:  I  /32*  NLQM  system  with  assimilation  of  altimeter  track 
data  from  ERS*2>  GFO  and  Jason- 1  altimeters. 

1/16°  NLOM:  Operational  1/16°  NLOM  system  with  assimilation  of 
real-time  altimeter  data  from  the  ERS-2  and  GFO  altimeters  (Jason- 1 
not  in  the  operational  data  stream  at  that  time). 

1/32°  NLOMn:  1/32°  NLOM  with  no  assimilation  of  ocean  data*  only 
atmospheric  forcing. 

NP:  eddy  not  present. 

*A  single  “fused"  model  eddy  represents  two  eddies  in  the  ocean  color 
and  the  model  eddy  center  lies  between  the  two  observed  eddies.  The 
model  eddy  center  position  error  is  listed  only  under  the  cIosct  ocean 
color  eddy  in  Table  I ,  Eddies  16  and  19  are  fused  in  both  1/32°  NLOM 
and  1/32°  NLOMn,  with  eddy  center  location  emirs  of  23  km  (26  km) 
for  eddy  #19  and  53  km  (54  km)  for  eddy  #16  in  1/32°  NLOM  <1/32* 
NLOMn).  Eddies  11  and  1 5  are  fused  in  1/16*  NLOM  with  eddy  center 
location  errors  of  55  km  for  eddy  #15  and  60  km  for  eddy  #11.  See 
Section  3.L4  for  discussion  of  the  fused  eddies. 

Both  1/32*  model  versions  of  eddy  #18  lie  east  of  63 °E  and  thus  are 
not  depicted  in  Fig,  2B,  D  or  Fig,  3B, 


forcing  is  sufficient  to  generate  eddy  #  14,  Fig*  3B,  D  shows 
that  a  Jason- 1  track  (not  assimilated  by  1/16°  NLOM)  also 
influenced  the  depiction  of  eddy  #14  in  1/32°  NLOM*  but 
otherwise  eddy  #14  is  poorly  observed  by  the  most  recent 
3  w  eeks  of  satellite  altimetry. 

Overall  1/32°  NLOM  captured  80%  of  the  10 
smallest  eddies  in  Table  1  vs  60%  by  1/16°  NLOM  and 
50%  by  1/32°  NLOMn  without  the  assimilation,  and  1/ 
32°  NLOM  had  the  lowest  eddy  center  position  error  for 
55%  of  these  eddies  vs  25%  for  1/16°  NLOM  and  20% 
for  1/32°  NLOMn,  Further,  the  median  position  error,  for 
the  portion  of  the  10  smallest  eddies  depicted,  is  22,5  km 
for  the  1/32°  NLOM,  32,5  km  for  the  1/16°  NLOM  and 
42  km  for  1/32°  NLOMn.  Among  these  eddies,  the 
smallest  location  error  is  1 1  km  for  1/32°  NLOM,  14  km 
for  1/16°  NLOM  and  26  km  for  1/32°  NLOMn.  Five  of 
the  eight  small  eddies  depicted  by  1/32°  NLOM  and  two 
of  the  six  in  1/16°  NLOM  have  location  errors  <  26  km 
and  87.5%  of  the  small  eddies  in  1/32°  NLOM  and  83% 
of  those  in  1/16°  NLOM  have  lower  eddy  center  position 
error  than  the  1/32°  NLOMn  median  of  42  km, 

3.1.3.  Source  of  improvement:  model  resolution  or 
number  of  altimeters? 

The  preceding  comparisons  between  the  data-assim- 
ilative  and  non-ass  imitative  models  clearly  demonstrate 
the  effectiveness  of  altimeter  SSH  assimilation  in  map¬ 
ping  ocean  eddies,  including  the  smaller  eddies*  In  ad¬ 
dition,  the  1/32°  NLOM  with  the  assimilation  clearly 
outperforms  the  1/16°  NLOM  with  assimilation.  Howev¬ 
er,  without  additional  analysis  it  is  not  clear  whether  this  is 
due  primarily  to  the  increased  model  resolution  or  to  the 
increase  from  2  to  3  altimeters  used  for  SSH  assimilation. 
We  do  know  that  the  1/32°  model  can  depict  eddies  with 
diameters  >  1/4°  vs  >  1/2°  for  the  1/16°  model.  All  four  of 
the  “1  /4°”  eddies  ( 1 7-20)  are  depicted  in  the  1  /32°  model 
(three  with  accurate  eddy  center  locations),  while  the  1/ 
1 6°  model  captures  only  one  of  these  eddies  (by  depicting 
it  as  an  eddy  >  1/2°  in  diameter). 

This  issue  is  investigated  further  using  Fig.  3B,  C,  D 
with  the  altimeter  tracks  from  all  3  altimeters  overlaid  on 
the  1/32°  and  1/16°  dataassimilative  results.  The  Jason-1 
altimeter  tracks  (in  white)  are  included  for  reference 
even  though  the  1/16°  model  did  not  assimilate  them 
(Fig.  3C).  In  all,  6  plots  like  Fig*  3B-D  were  examined, 
22  Sept,,  29  Sept,  and  6  Oct.  from  the  data-assimilative 
1/16°  and  1/32°  NLOM.  Each  had  the  most  recent  week 
of  altimeter  tracks  overlaid  as  solid  lines  (the  rest  as 
dashed  lines).  In  total,  these  cover  the  most  recent 
3  weeks  of  altimeter  data  (16  Sept*-6  Oct.)*  Although  a 
week  of  altimeter  tracks  were  overlain  as  solid  lines  on 
each  plot,  the  analyses  were  performed  daily  using  a  3- 
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day  data  window.  The  length  of  time  used  for  examining 
earlier  results  is  based  on  the  model  time-scale  for 
forecast  skill  in  this  region,  as  shown  in  Section  3.5*  The 
influence  of  the  earlier  results  on  the  results  and  com¬ 
parisons  for  6  Oct.  were  dearly  evident. 

Examination  of  the  plots  indicated  that  for  80%  of  the 
20  eddies,  the  existence,  location  and  amplitude  in  the 
models  were  not  noticeably  influenced  by  this  particular 
configuration  of  observed  Jason- 1  tracks.  Eddies  9  and 
14  were  clearly  influenced  and  eddies  6  and  8  may  have 
been  influenced.  Despite  the  Jason- 1  influence  on  eddy 
#14  in  1/32°  NLOM,  the  results  were  not  improved,  as 
discussed  earlier  and  shown  in  Table  L  In  the  case  of 
eddy  #9,  a  recent  Jason- 1  track  crossed  the  interior  of  the 
cyclonic  eddy  depicted  by  the  1/32°  NLOM  on  6  Oct. 
(Fig.  3B),  a  location  where  the  model  shows  a  SSH  ridge 
between  two  cyclonic  eddies  on  29  Sept,  (Fig,  3D),  In 
contrast,  an  anti-cyclonic  eddy  was  present  in  1/16° 
NLOM,  which  did  not  assimilate  the  Jason- 1  data 
(Fig.  3C).  Because  the  29  Sept,  cyclonic  eddy  east  of 
the  SSH  ridge  in  1/32°  NLOM  was  propagating  west¬ 
ward,  it  is  quite  likely  it  would  have  merged  with  the 
eddy  west  of  the  ridge  even  without  assimilation  of  the 
Jason- 1  track.  However,  that  track  would  have  clearly 
benefited  1/16°  NLOM.  This  raises  the  issue  of  dyna¬ 
mical  interpolation  skill  and  the  impact  of  the  data 
assimilation  earlier  in  time.  At  the  beginning  of  the  time 
frame  examined  (2  weeks  before  6  Oct),  the  1/32° 
NLOM  had  already  pre-conditioned  itself  with  cyclonic 
eddies  in  the  Gulf  of  Oman,  while  the  1/16°  NLOM  was 
already  pre-conditioned  with  anti -cyclonic  eddies.  The 
1/32°  NLOM  without  the  data  assimilation  shows  only 
very  weak  features  in  the  Gulf  of  Oman  (Fig.  2D), 

The  preceding  results  strongly  indicate  that  the  in¬ 
crease  in  model  resolution  had  a  larger  positive  impact  on 
the  results  than  the  addition  of  the  Jason- 1  tracks  in  this 
particular  situation.  In  general,  earlier  studies  have  found 
that  adding  a  third  nadir  beam  altimeter  modestly  reduces 
the  RMS  SSH  analysis  error  (Le  Traon  and  Dibarboure, 
1999;  Hurlburt  et  al.,  2000;  Le  Traon  et  aL,  2001; 
Smedstad  et  aL*  2003).  Outside  of  Section  3,1 ,  SSH  from 
the  same  number  of  altimeters  was  assimilated  by  1/16° 
and  1/32°  NLOM  in  most  of  the  tests  and  all  of  the 
comparisons  of  real-time  results.  In  Section  3,4,  the 
Hawaiian  Islands  region  is  used  to  compare  results  where 
both  models  assimilated  data  from  one  or  three  altimeters. 

3.1.4,  The  roles  of  altimeter  data  assimilation  and 
ocean  model  dynamical  interpolation/forecast  skill  in 
nowcasting  small  eddies  at  1/16°  and  1/32°  resolution 
Three  noteworthy  examples  illustrate  the  importance 
of  ocean  model  dynamical  in terp-olaii on/forecast  skill 


and  the  value  of  increased  model  resolution  in  now- 
casting  small  eddies.  The  pair  of  counter-rotating  eddies, 
#12  and  #13,  provide  an  excellent  example.  Also  note  the 
current  that  enters  the  region  near  6Q°E  and  ultimately 
wraps  around  eddy  #1 2.  It  is  marked  by  a  ribbon  of  blue 
in  the  ocean  color  and  flows  northward  along  -  60QE, 
then  westward  north  of  18°N  and  then  northward  along 
-  58.5°E  until  it  wraps  around  eddy  #12,  The  two  eddies 
and  the  current  are  seen  in  both  the  1/16°  and  1/32° 
NLOM  with  assimilation  and  they  are  depicted  with 
similar  accuracy  in  comparison  to  the  ocean  color,  except 
that  the  center  location  of  eddy  #13  is  more  accurate  in 
1/32°  NLOM.  Neither  of  the  two  eddies  nor  the  current  is 
seen  in  1/32°  NLOMn  with  no  assimilation.  However,  it 
is  evident  in  Fig,  3B,  C  that  eddies  #12  and  #13  were 
unobserved  by  the  most  recent  week  of  altimeter  data  and 
only  the  southernmost  portion  of  the  current  was  ob¬ 
served.  Hence,  we  must  explore  earlier  in  time  to  in¬ 
vestigate  the  role  of  the  altimeter  data  in  forming  these 
observed  features  and  the  dynamical  interpolation  and 
forecast  skill  of  the  models  in  creating  a  realistic  de¬ 
piction  on  6  Oct. 

To  start,  we  note  that  within  the  pictured  region,  the 
length  of  the  current  feeding  into  eddy  #12  is  -  600  km. 
Two  and  one  half  weeks  would  be  required  for  a  40  cm 
s  1  current,  comparable  to  the  modelled  current  speed,  to 
advect  the  low  chlorophyll  water  that  distance,  implying 
that  at  least  the  southern  part  of  the  current  has  been  a 
persistent  feature  for  at  least  that  long.  Going  back 
1  week  (altimeter  data  between  I  and  2  weeks  prior  to 
6  Oct.  overlaid),  we  see  that  the  current  was  better  ob¬ 
served  but  the  two  eddies  were  only  peripherally  ob¬ 
served  (Fig,  3D),  Tlie  current  and  eddy  #1 2  were  already 
in  place.  Eddy  #!  3  was  just  forming  in  1/32°  NLOM  but 
was  already  in  place  in  \/\6°  NLOM  (not  shown).  Going 
back  2  weeks  (altimeter  data  2  to  3  weeks  prior  to  6  Oct. 
overlaid)  (not  shown),  we  find  the  genesis  of  eddy  #  1 3  as 
a  WNW  bulge  from  eddy  #2.  This  bulge  is  well  observed 
by  2  ERS-2  and  a  GFO  track  and  peripherally  by  a  Jason- 
1  track,  Eddy  #13  is  not  yet  present  in  1/32°  NLOM  and 
is  just  forming  in  1/16°  NLOM.  Thus,  altimeter  data 
2-3  weeks  prior  was  crucial  to  the  depiction  of  eddy  #  1 3 
on  6  Oct,  possibly  with  some  assistance  in  the  eddy  pinch 
off  from  a  GFO  track  a  week  later  (Fig.  3D).  Eddy  #12 
was  already  present  on  22  Sept,  but  only  peripherally 
observed.  Thus,  this  eddy  persisted  for  at  least  2  weeks 
with  little  change  and  was  only  peripherally  observed  by 
3  weeks  of  altimeter  data,  but  still  it  is  realistically 
represented  in  1/16°  and  1/32°  NLOM  (but  not  in  1/32° 
NLOMn  with  no  assimilation).  A  shifted  version  of  the 
current  was  present  on  22  Sept,  and  it  was  better 
observed  than  on  29  Sept,  or  6  Oct. 
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The  second  example  is  the  duster  of  6  small  eddies 
consisting  of#  19,  16,  IS,  17,  11  and  15  south  to  north, 
where  #18  is  the  only  ami-cyclonic  eddy  in  this  group. 
The  entire  cluster  spans  just  over  2°N-3  and  1°E~W.  All 
of  these  eddies  are  captured  by  at  least  one  of  the  two  data 
assimilative  systems,  4  by  1/32°  NLOM  and  3  by  1/16° 
NLOM,  but  1/32°  NLOMn  with  no  data  assimilation 
also  captured  4  of  the  eddies  (Table  1)*  However,  for  the 
eddies  in  the  cluster,  the  eddy  center  location  error  is  less 
than  the  smallest  location  error  in  1/32°  NLOMn  (26  km) 
for  75%  of  the  1/32°  and  67%  of  the  1/!  6°  NLOM  eddies 
depicted.  This  is  consistent  with  the  overall  results  for  the 
1 0  smallest  eddies  presented  earlier,  which  showed  value 
added  for  altimeter  data  assimilation  in  representing 
small  eddies  25-75  km  in  diameter.  By  examining  the 
six  plots  listed  earlier,  it  was  evident  that  the  model 
depictions  of  all  the  eddies  were  influenced  by  the 
altimeter  data,  mainly  from  ERS-2  and  GFO.  The  duster 
was  only  minimally  and  peripherally  observed  by 
Jason- L 

Starting  with  the  southernmost  eddy  and  progressing 
northward,  eddy  #19  was  well  observed  by  an  ERS-2 
track  during  the  week  leading  up  to  22  Sept.  The  result 
on  22  Sept,  was  a  small  strong  eddy  in  1/32°  NLOM  and 
in  1/16°  NLOM  a  weaker  eddy  and  narrow  SSH  trough 
connected  to  eddy  #L  During  the  next  2  weeks  (through 
6  Oct),  eddy  #19  was  not  observed  by  the  altimetry.  It 
was  well  maintained  by  1/32°  NLOM,  but  the  weak  eddy 
and  sharp  trough  in  1/16°  NLOM  disappeared  by  29 
Sept.  During  the  week  leading  up  to  29  Sept.,  the 
location  of  nearby  eddy  #  1 6  on  6  Oct.  was  well  observed 
at  the  crossover  of  ERS-2  and  GFO  tracks,  but  no  eddy 
resulted  in  either  model.  In  1/32°  NLOM  there  was  an 
ENE  extension  of  the  SSH  trough  associated  with  eddy 
#19  which  had  the  ERS-2  track  crossing  the  eastern  edge 
of  it  just  north  of  the  crossover.  During  the  week  leading 
up  to  6  Oct,  a  GFO  track  crossed  the  center  of  eddy  #  1 6 
and  an  ERS-2  track  crossed  the  eastern  interior.  The 
result  was  a  weak  eddy  #1 6  in  the  1/16°  NLOM  (14-km 
eddy  center  location  error),  but  only  a  modification  of  the 
existing  SSH  trough  in  1/32°  NLOM.  In  1/32°  NLOM 
eddies  16  and  19  were  fused  with  the  model  eddy 
encompassing  both  and  with  the  eddy  center  between 
them,  but  much  closer  to  eddy  #19  (23-km  eddy  center 
location  error).  Anti-cyclonic  eddy  #18  just  to  the  north 
was  well  observed  by  altimetry  during  the  weeks  leading 
up  to  both  29  Sept,  and  6  Oct.  Both  1/16°  and  1/32° 
NLOM  were  pre-conditioned  to  represent  this  eddy,  but 
1/16°  NLOM  was  better  pre-conditioned  and  matched 
the  observed  eddy  center  much  more  accurately. 

Northwest  of  eddy  #  1 8  an  ERS-2  track  passed  close  to 
the  observed  location  of  eddy  #17  on  6  Oct.  during  the 


week  leading  up  to  29  Sept,  This  formed  an  SSH  trough 
linked  to  eddy  #1  in  1/32°  NLOM,  During  the  week 
leading  up  to  6  Oct.  a  GFO  track  crossed  the  neck  of  the 
trough  and  also  formed  an  SSH  ridge  to  the  south 
narrowing  the  trough,  both  assisting  in  cutting  off  eddy 
#17  and  giving  an  eddy  center  location  error  of  only 
1 3  km.  The  1/16°  NLOM  did  not  represent  these  features 
with  sufficient  amplitude  and  sharpness  to  form  eddy 
#17,  which  had  an  SSH  diameter  of  only  1/2°E-W  by 
1/4°N-S  in  1/32°  NLOM,  smaller  than  the  1/2°  limit  for 
realistic  representation  of  eddies  in  1/16°  NLOM, 

An  ERS-2  and  a  GFO  track  passed  close  to  the  enter 
of  eddy  #  1 5  during  the  week  leading  up  to  6  Oct,  but  did 
not  create  or  pass  over  an  eddy  in  either  model  (#15  was 
not  present  in  1/32°  NLOM).  A  GFO  track  crossed  the 
eddy  in  1/16°  NLOM  during  the  week  leading  up  to 
22  Sept,  and  the  eddy  propagated  westward  toward  the 
observed  location  of  eddy  #15  on  6  Oct,  The  1/32° 
NLOM  was  not  well  pre-conditioned  to  represent  this 
eddy.  Slightly  farther  to  the  south  (and  east  of  63°E),  the 
same  GFO  track  crossed  the  northeastern  quadrant 
(southwestern  quadrant)  of  a  larger  cyclonic  eddy 
depicted  in  1/32°  NLOM  (1/16°  NLOM)  on  22  Sept. 
During  the  week  leading  up  to  29  Sept.,  a  GFO  track 
crossed  the  southeastern  quadrant  of  the  eddy  in  1/32° 
NLOM,  favorably  impacting  its  westward  propagation 
into  the  observed  position  of  eddy  #11  on  6  Oct.  (22  -km 
eddy  center  location  error).  This  occurred  during  a  week 
when  it  was  not  observed  by  altimetry.  The  eddy  in  l/l  6s 
NLOM  lay  east  of  the  GFO  track  on  29  Sept.  Thus,  its 
propagation  did  not  benefit  from  the  GFO  track  and  eddy 
#11  was  not  represented  in  1/1 6°  NLOM, 

The  smallest  eddy,  #20,  provides  the  third  example.  It 
is  cyclonic  and  was  not  observed  by  altimetry.  It  formed 
where  a  boundary  current  associated  with  anti-cyclonic 
eddy  #3  separated  from  the  boundary,  Eddy  #3  and  the 
boundary  current  separation  are  represented  with  com¬ 
parable  realism  in  both  1/16°  NLOM  and  1/32°  NLOM. 
However,  eddy  #20  is  seen  only  in  1/32°  NLOM.  At  1/4° 
in  diameter  in  both  the  model  and  the  ocean  color,  it  is 
smaller  than  the  1/2°  minimum  for  realistic  depiction  in 
1/16°  NLOM,  In  1/32°  NLOM  three  such  eddies  are 
observed  farther  south  on  22  Sept,  (not  shown).  As  eddy 
#3  moves  northward,  1/32°  NLOM  depicts  a  cyclonic 
eddy  (#20)  in  the  location  observed  in  the  ocean  color  on 
6  Oct,  (Fig.  3A)  on  both  29  Sept.  (Fig.  3D)  and  6  Oct. 
(Fig.  3B)  (12-km  eddy  center  location  error  on  6  Oct.). 
The  accurate  location  of  the  eddy  may  be  facilitated  by  a 
slight  boundary  indentation  not  present  in  the  1/16° 
NLOM  boundary. 

The  preceding  results  show  that  1/16°  and  1/32° 
NLOM  are  skilled  at  using  multiple  altimeter  tracks  from 
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multiple  satellites,  assimilated  through  time  and  space,  to 
build  coherent  mesoscale  features  that  are  consistent 
with  observations,  including  small  eddies.  They  also 
demonstrate  that  1/32°  NLGM  has  greater  skill  than 
1/16°  NLOM,  especially  for  small  eddies  25-75  km  in 
diameter.  In  addition,  they  show  that  a  model  with  dy¬ 
namical  interpolation  skill  can  use  data  assimilation  to 
produce  an  ambience  more  conducive  to  the  formation  of 
poorly  observed  eddies  than  the  same  model  without  the 
assimilation,  a  capability  that  is  especially  useful  in  rep¬ 
resenting  small  eddies  and  one  that  is  clearly  demonstrat¬ 
ed  by  the  results  in  this  subsection. 

3.2  Caribbean/Gulf  of  Mexico 

The  increase  in  horizontal  resolution  produces  sig¬ 
nificant  changes  in  the  Caribbean/Gulf  of  Mexico  region 
{Fig.  4).  Transport  through  II  of  the  14  Caribbean 
passages  is  more  realistic  at  the  finer  resolution,  in¬ 
cluding  more  realistic  Gulf  of  Mexico  inflow  through  the 


Yucatan  Channel.  The  improvement  in  the  Yucatan 
Channel  transport  occurs  primarily  because  the  flow  in 
the  1/32°  model  is  more  surface-trapped  than  at  1/16°,  so 
more  of  the  transport  passes  through  the  Caribbean 
passages  {since  more  is  above  the  sill  depths).  In  the  1/ 
32°  model  there  is  less  transport  through  the  Northwest 
Providence  Channel  (feeds  into  the  western  boundary 
current  near  26DN)  in  better  agreement  with  observa¬ 
tions.  Unlike  the  1/16°  model,  the  narrowest  point  in  the 
Florida  Strait  (along  80°W  extending  back  to  the  Gulf  of 
Mexico  outflow)  is  north  of  this  channel,  as  observed. 

The  improvement  in  Yucatan  Channel  transport  is 
particularly  significant  because  it  feeds  the  Loop  Current 
which  sheds  large  eddies  in  the  Gulf  of  Mexico.  The 
value  added  of  the  finer  resolution  in  the  Gulf  of  Mexico 
is  illustrated  in  a  comparison  of  SeaWiFS  ocean  color 
with  SSH  from  the  1/16°  and  1/32°  models  (Pig.  5),  Both 
assimilate  SSH  data  from  the  GFO  and  Jason- 1  alti¬ 
meters  plus  SSI.  The  key  features  are  the  Loop  Current 
and  a  shed  eddy.  The  SeaWiFS  imagery  shows  these  as 
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Fig.  4.  Observed  transports  in  Sv  ( 10*  m3  s  ’)  vs  model  mean  (1993-2000)  transports  through  Intra- Americas  Sea  passages  at  resolutions  of  1/16° 
(with  data  assimilation)  and  1/32*  (without  data  assimilation),  the  most  suitable  estimates  we  had  from  each  model  over  an  extended  common  lime 
period.  There  arc  no  data-assimilative  1/32°  model  results  prior  to  Oct.  2000  and  the  passage  transports  are  not  strongly  constrained  by  the  data 
assimilation.  Also  shown  are  6-year  mean  currents  from  the  1/32°  model  (climatological  forcing  only)  superimposed  on  their  strength.  Johns  ct  al. 
(2002)  provided  the  best  estimates  based  on  observations  and  did  not  adjust  them  to  provide  a  closed  mass  budget.  The  “observed"  value  of  29.2  Sv 
for  the  Yucatan  Channel  used  here  was  obtained  by  subtracting  the  two  observed  Florida  Strait  inflows  between  Yucatan  and  27°N  (1.2  and  1 ,9  Sv) 
from  the  well  determined  value  of  32.3  Sv  at  27*N  given  by  Larsen  (1992).  The  sum  from  the  remaining  Caribbean  passages  gives  a  second 
observationally  based  estimate  of  28.0  Sv  for  Yucatan  Channel.  [Figure  courtesy  T.L.  Townsend,  Naval  Research  Laboratory]. 
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Fig.  5.  A  comparison  of  SoaWiFS  ocean  color  with  SSH  from  the  1/32°  and  1/16°  data  assimilative  systems.  Model  SSH  from  28  Aug.  2003  overlain 
on  SeaWiFS  chlorophyll  concentration  from  22  to  28  Aug.  2003  latest  cloud  free  pixel  composite.  SSH  data  from  the  GFO  and  Jason- 1  altimeters 
were  assimilated,  SSH  contours  >50  cm  arc  omitted. 


1/16°  Global  NLOM 


dark  areas  of  low  chlorophyll  outlined  by  a  Mississippi 
River  plume  of  high  chlorophyll  (red  and  yellow)  along 
the  eastern  edge  of  the  eddy  and  the  Loop  Current.  The  1  / 
16°  model  largely  misses  the  shed  eddy.  In  an  animation 
covering  June-Sept  2003,  the  1/32°  NLOM  routinely 
outperformed  the  1/16°  model.  In  a  more  detailed  com¬ 
parison  of  5  data-assimilative  ocean  prediction  systems 
covering  the  Gulf  of  Mexico,  Chassignet  et  ah  (2005) 
found  that  the  1/3 2a  NLOM  system  performed  the  best 
in  comparison  to  ocean  color  imagery  with  a  1/12° 
Atlantic  HYCOM  system  a  close  second.  Animations 
comparing  these  five  systems  can  be  found  on  the 
NLOM  and  HYCOM  web  pages  (http: //www. ocean. 
nrlssc.navy.mil/giobaLnlom  and  http://www.hy com. 
org).  Besides  the  improved  inflow  through  Yucatan 
Channel,  other  reasons  for  better  performance  at  1/32° 
resolution  are  better  dynamical  interpolation  of  the 
assimilated  altimeter  data  and  elimination  of  a  southward 
leakage  of  the  Deep  Western  Boundary  Current  through 
the  Florida  Straits  (found  in  the  1/16°  model)  via  a 
modification  to  the  bottom  boundary  condition.  The 
southward  leakage  in  the  1/16°  model  resulted  in  a 
westward  subthennocline  current  along  the  continental 
slope  in  the  northern  Gulf  of  Mexico.  Such  a  current  can 
advect  upper  ocean  currents  (Hurlburt  and  Thompson, 


1980;  Hurlburt  et  al.,  1996)  (e.g.  the  semi-detached  Loop 
Current  eddy  in  Fig.  5),  and  in  this  case  damage  the 
dynamical  interpolation  skill  of  the  1/16°  model  in  the 
northern  Gulf  of  Mexico,  as  illustrated  in  the  figure. 

33.  Gulf  Stream 

Fig.  6  shows  results  from  the  c 1  i mato logically- forced 
1/32°  global  model  in  the  Gulf  Stream  region,  one  of  the 
most  challenging  regions  to  simulate  in  the  world  ocean. 
It  also  shows  comparisons  to  the  corresponding  1/1 6° 
climatological  global  simulation  and  to  observations. 
The  mean  SSH  from  each  model  is  overlain  by  the 
\  5-year  mean  north  wall  of  the  Gulf  Stream  from  satellite 
IR  (middle  line)  flanked  by  lines  one  standard  deviation 
away.  The  RMS  SSH  variability  is  compared  with  the 
Gulf  Stream  mean  axis  (center)  flanked  by  the  standard 
deviation  about  the  mean  axis  and  extreme  positions 
determined  from  Topex/Poseidon  data  (Lee,  1997).  Both 
models  show  a  realistic  mean  Gulf  Stream  pathway, 
although  the  1/16°  model  has  a  slight  northward  bulge 
near  70°W.  Resolution  of  1/16°  or  —  7  km  is  the  min¬ 
imum  required  to  obtain  a  realistic  Gulf  Stream  pathway 
(Hurlburt  and  Hogan,  2000),  At  coarser  resolution,  the 
pathway  is  very  unrealistic  and  the  Gulf  Stream  transport 
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Fig.  6.  Mean  SSH  (A,  C)  and  RMS  SSH  variability  (B,  D)  in  the  Gulf  Stream  region  from  the  1/16°  and  1/32°  global  NLOM  with  climatological 
forcing  and  no  assimilation  of  SSTor  SSH  data.  Panel  E  depicts  RMS  SSH  variability  ( 1993-2003)  computed  from  the  model  independent  MODAS 
2D  daily  altimetric  SSH  analyses.  The  lines  overlain  in  A  and  C  depict  the  mean  pathway  of  the  Gulf  Stream  1R  north  wall  L  9B2—  1  996  ±  1  standard 
deviation  by  Conation  and  Sirfces  (unpublished  work).  The  5  lines  overlain  in  B,  D  and  E  are  the  Gulf  Stream  mean  axis  (center),  standard  deviation 
about  the  mean  axis  and  extreme  positions  determined  from  Topex/Poseidon  data  (Lee,  1997).  The  Gulf  Stream  (depicted  by  the  lightly  packed 
contours)  robustly  separates  from  the  coast  near  Cape  Hatteras  (approximately  363N,  75*  W)  and  follows  a  realistic  pathway  to  a  region  south  of  the 
Grand  Banks  (approximately  42 DN,  50° W). 


is  much  too  weak.  However,  a  realistic  pathway  is  more 
easily  and  robustly  obtained  at  1/32°  resolution  and  the 
Gulf  Stream  transport  is  much  higher  and  more  realistic, 
particularly  south  of  the  Grand  Banks  at  50°W  between 
36  and  44°N  (77  Sv  at  1/32*  and  42  Sv  at  1/16°  for  the 
baroclinic  transport  vs  an  estimate  of  65  Sv  in  the  upper 
1000  m  by  Dietrich  et  ah  (1980)).  Associated  with  this 
large  increase  in  transport  is  a  large  increase  in  the 
strength  and  eastward  penetration  of  the  non-linear 
recirculation  gyre  on  the  south  side  of  the  Gulf  Stream 
(shown  in  red  and  black).  This  alters  the  large-scale  C- 
shape  of  the  subtropical  gyre. 


The  RMS  SSH  variability  is  more  realistically  dis¬ 
tributed  in  the  1/32°  model,  with  higher  variability  pro¬ 
truding  farther  to  the  east,  consistent  with  RMS  SSH 
(1993-2003)  computed  from  the  model-independent 
MODAS2D  SSH  analyses  that  are  operational  at 
NAVOCEANO  (Fox  et  al,  2002).  Further,  the  pattern 
and  amplitude  of  the  variability  are  helpful  in  assessing 
the  ability  of  the  climatologically-forced  model  (without 
data  assimilation)  to  simulate  specific  features  of  the 
mean  flow  and  variability.  Both  the  1/16°  and  1/32°  models 
show  the  narrow  corridor  of  Gulf  Stream  variability  west  of 
-  69°W  with  low  variability  to  the  north  of  this  corridor 
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However,  the  corridor  bulges  too  far  north  m  the  1/1 6° 
model.  Both  models  also  simulate  the  southward  ex¬ 
tension  of  high  variability  near  the  New  England  Sea¬ 
mount  Chain  60°-65°W),  a  region  of  high  cold  core 
ring  generation  with  a  westward  trail  of  elevated  vari¬ 
ability  as  the  rings  propagate  to  the  west,  features  better 
represented  in  the  1/32°  model.  The  largest  advantage  of 
the  1/32°  model  is  demonstrated  by  using  the  pattern  of 
high  variabi  lity  to  assess  the  eastward  penetration  of  the 
non-linear  recirculation  gyre  and  its  latitudinal  position 
in  that  region.  Both  the  1/32°  model  and  MGDAS2D 
variability  show  eastward  penetration  to  ~  44°W  near 
37°-38°N  and  a  wraparound  pattern  of  high  variability 
to  the  southern  edge  of  the  gyre  along  35°-36°R  How¬ 
ever,  even  the  1/32°  model  fails  to  reproduce  the  north¬ 
eastward  extension  of  high  variability  (41°-45°N,  40°- 
45°W)  associated  with  the  Mann  Eddy  east  of  the  Grand 


Banks.  See  Hurlburt  and  Hogan  (2000)  for  additional 
discussion  of  other  NLOM  simulations  vs  observed  SSH 
variability  in  the  Gulf  Stream  region. 

The  improvements  in  the  simulation  of  the  Gulf 
Stream  at  1/32°  resolution  have  nowcast/forecast 
implications,  as  is  illustrated  in  Fig.  7.  The  Gulf  Stream 
in  the  data  assimilative  1/32°  model  is  more  realistic, 
following  the  meanders  indicated  in  the  operational  JR 
frontal  analyses  performed  at  NAVOCEANO  better  than 
at  1/16°.  In  addition,  the  Gulf  Stream  is  stronger  and 
more  inertial  in  the  1/32°  model  as  indicated  by  the  larger 
change  in  SSH  across  the  stream.  The  improvements  in 
the  nowcast  Gulf  Stream  pathway  provide  a  better 
forecast  initial  condition,  and  hence  a  more  accurate 
forecast  as  well  This  is  seen  in  axis  error  statistics  shown 
in  Fig.  8.  The  error  is  calculated  as  the  area  between 
corresponding  SSH  contours  in  the  forecast  and  the 
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Fig.  7.  Nowcast  SSH  in  the  Gulf  Stream  region  from  the  (A)  and  (B)  1/16°  global  NLOM  systems  for  12  May  2003.  The  lines  overlain  are 
operational  1R  frontal  analyses  by  the  Naval  Oceanographic  Office;  black  lines  are  segments  based  on  IR  data  >4  days  old. 
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Forecast  length  (days) 


Fig.  8.  Mean  axis  error  for  24  Gulf  Stream  forecasts  from  the  1/16° 
(dashed  tines)  and  1/32°  (solid  lines)  nowcast/forccast  systems 
initialized  during  the  time  period  \  June  2001-31  May  2002.  The 
forecasts  from  the  1  /32°  model  have  lower  axis  error  for  the  duration  of 
the  forecast 

verifying  analysis  divided  by  the  contour  length  in  the 
analysis.  Since  we  are  using  the  model  analysis  as  the 
“truth”,  we  use  several  different  SSH  contours  to 
represent  the  Gulf  Stream  pathway  and  calculate  an 
average  error  for  these  contours.  Throughout  the  forecast 
period,  the  l  /32°  system  exhibits  lower  error  for  the  Gulf 
Stream  axis. 

3.4 ,  Hawaii 

The  ability  of  satellite  altimetry  to  constrain  an  eddy¬ 
resolving  model  was  discussed  in  Section  3.T  Where 
the  mesoscale  variability  is  highly  non-deterministic,  as 
in  the  region  around  the  Hawaiian  Islands,  insufficient 
data  to  constrain  the  model  can  introduce  significant 
errors  into  the  nowcast.  Hurlburt  et  ai.  (2000)  and 
Smedstad  et  al.  (2003)  discuss  the  suppression  of 
nowcast  SSH  error  as  the  number  of  altimeters  used  in 
assimilation  is  increased  Fig.  9  (A)  and  (B)  shows  a 
comparison  of  1/16°  and  1/32°  NLOM  SSH  snapshots 
for  19  March  2004  in  the  Hawaii  region  when  data  from 
only  1  satellite  altimeter  (Jason- 1 )  were  assimilated.  The 
snapshots  show  some  qualitative  similarities,  but  large 
differences  can  also  be  noted  (e.g.,  robust  eddies  in  the 
1/32°  SSH  nowcast  at  approximately  157°W/18.5°N, 
157°W/19.8°N  and  l56°W/22dN  that  the  I/I6&  model 
doesn't  exhibit).  This  occurs  because  data  from  1 
altimeter  are  not  sufficient  to  adequately  constrain  the 
models  in  this  region  and  thus  they  are  unable  to 
skillfully  map  the  mesoscale  eddies,  which  are  predom¬ 
inantly  non-deterministic  in  relation  to  atmospheric 
forcing. 


Fig.  9  (C)  and  (D)  show  a  comparison  of  SSH 
snapshots  from  the  1/32°  and  1/16°  models  on  5  June 
2004  when  data  from  3  satellites  (Jason- U  GFO  and 
Envisat)  were  assimilated.  The  solutions  show  better 
agreement  with  an  anomaly  correlation  of  .94  for  the 
SSH  in  panels  (C)  and  (D)  vs  .30  for  panels  (A)  and  (B). 
In  addition,  animations  of  this  area  from  the  1/32°  model 
tend  to  show  better  time  continuity  of  small-scale 
features.  In  the  Hawaii  region  we  see  much  larger 
impacts  from  the  resolution  increase  and  the  data 
increase  from  one  to  three  altimeters  than  in  the  northern 
Arabian  Sea  (Figs.  2,  3,  Table  1).  In  that  case  the 
resolution  increase  and  an  increase  from  2  to  3  altimeters 
gave  clear  improvement,  but  less  dramatic  than  here. 

The  finer  resolution  model  also  gives  improvement 
in  comparison  to  SSH  from  four  tide  gauges  in  the 
Hawaiian  Islands  (30-day  moving  average  applied,  not 
shown)  with  reduced  RMS  error  at  all  four  and  an 
overall  19%  reduction  in  the  error  variance  during  the 
1  -year  period  (June,  200 1  -  May,  2002)  used  for  compar¬ 
ison  and  statistical  evaluation  of  the  two  systems. 
During  this  period  data  from  3  altimeters  (ERS-2,  GFO 
and  Topex/Posetdon)  were  consistently  available.  It  is 
noteworthy  that  while  the  error  variance  was  reduced  by 
1 9%,  the  SSH  variance  at  all  4  tide  gauge  locations  (over 
the  same  time  period)  increased  in  comparison  to  the  1/ 
16°  model,  overall  by  2.2*  when  the  1/32°  model  was 
used  to  assimilate  the  data.  This  improves  agreement 
with  the  tide  gauge  variance,  i.e.  the  variance  is  4.3  cm2 
for  the  1/16°  model,  9.7  cm2  for  the  1/32°  model,  and 
1 5  cm2  for  the  tide  gauge  data. 

5.5.  Evaluation  of  nowcast/forecast  skill 

5.5. /.  Model-data  comparisons  using  moored  buoys 
and  tide  gauges 

As  in  Smedstad  et  ah  (2003),  modeled  SSH  is 
compared  with  observations  using  daily  time  series  from 
94  tide  gauge  stations.  A  30-day  running  mean  was 
applied  to  filter  time-scales  not  well  observed  by 
altimetry.  Fig.  10  shows  a  histogram  of  RMS  difference 
and  correlation  statistics  between  model  nowcasts  and 
observed  sea  level  from  each  model.  It  is  important  to 
note  that  NLOM  assimilates  SSH  from  real-time 
satellite  altimetry,  not  tide  gauges.  Both  models  have 
comparable  RMS  difference  and  correlation  statistics 
(median  RMSD  ^4.0  cm  and  median  correlation  -  .80), 
The  distribution  of  correlations  is  nearly  identical,  but 
the  distribution  of  RMSD  shows  improvement,  with 
shifting  of  the  distribution  to  lower  RMSD  values. 

Fig.  11  shows  a  histogram  of  RMS  difference  and 
correlation  statistics  between  daily  time  series  of 
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Hawaiian  Islands  subregion 


Assimilation  of  SSH  data  from  1  altimeter,  19  Mar  2004 


Assimilation  of  SSH  data  from  3  altimeters,  S  Jun  2004 
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Fig.  9,  Snapshots  of  SSH  (with  surface  currents  overlain)  in  the  Hawaiian  islands  subregion  for  19  Mar  2004  (A  and  B)  and  5  Jun,  2004  (C  and  D) 
from  the  1/32°  (A,  C)  and  1/16*  (B,  D)  global  models  (both  with  assimilation).  For  19  Mar.  2004*  data  from  only  one  satellite  altimeter  (Jason- 1 )  were 
assimilated.  For  5  Jun.  2004,  data  from  3  satellite  altimeters  (Jason- 1 1  GFO  and  Envisat)  were  assimilated.  To  facilitate  inter-comparison,  a  constant 
offset  {4  cm)  was  subtracted  from  (B)  and  (D)  to  better  reference  them  to  the  areal  mean  in  (A)  and  (C), 


observed  SST  and  nowcast  SST  from  the  1/32°  and  1/ 
16°  models.  A  total  of  84  buoys  (from  the  National  Data 
Buoy  Center  (NDBC)  and  the  TAO  and  PIRATA  arrays) 
were  used  in  this  comparison.  NLOM  assimilates  the 
MODAS  daily  operational  SST  analyses  of  satellite 
MCSST  data  (Sections  2.2  and  2.5).  Both  models  have 
comparable  RMS  difference  and  correlation  statistics 
(median  RMSD«,35°C  and  median  correlation  ,94)* 
This  negligible  difference  is  due  to  the  strong  SST 
assimilation  used  in  both  models  and  lower  sensitivity  to 
horizontal  resolution  than  SSH. 


1.5.2  Forecast  skill 

Several  metrics  are  used  routinely  to  assess  model 
forecast  skill,  all  based  on  comparison  to  a  verifying 
analysis,  and  forecast  skill  superior  to  persistence  (a 
forecast  of  no  change),  and  either  climatology  or  some 
minimum  value  of  the  metric.  The  metrics  are  RMS 
error,  anomaly  correlation,  skill  score,  and  for  the 
Kuroshio  and  the  Gulf  Stream,  axis  error  (Fig.  8).  See 
Smedstad  et  al.  (2003)  for  discussion  and  examples  of 
results  from  each  one.  See  the  real-time  NLOM  web 
page  http://www.ocean.nrlssc.navy.mil/globaLnlom  for 
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A .  1/32°  NLOM  B.  1/16*  NLOM 
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Fig.  10,  Histogram  of  RMS  difference  (A  and  B}  and  correlation  statistics  (C  and  D)  between  observed  and  model  nowcast  sea  level  from  the  1/16° 
and  11/32*  NLOM  nowcast/forecast  systems  during  the  period  I  June  2001-31  May  2002.  Note  that  NLOM  assimilates  SSH  from  real-time  satellite 
altimetry,  not  tides  gauges.  These  comparisons  are  based  on  deviations  from  the  mean  of  each  time  series,  94  tide  gauge  stations  were  used. 
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Fig,  1 1 .  Histogram  of  RMS  difference  (A  and  B)  and  correlation  statistics  (C  and  D)  between  observed  SST  from  84  moored  buoys  and  nowcast  SST 
from  the  1/16°  and  1/32°  NLOM  now  cos  V  forecast  systems  during  the  period  1  June  2001-31  May  2002,  Note  that  NLOM  assimilates  SST  from  daily 
operational  MOD  AS  SST  analyses  of  satellite  MCSST  data,  which  do  not  include  data  from  moored  buoys. 
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real-time,  archived  and  summary  results  for  each  metric 
in  many  regions  from  both  the  1/16°  and  l//32°  models. 

Here  we  focus  on  one  metric,  anomaly  correlation 
(AC)  between  the  forecast  and  the  verifying  analysis,  as 
is  typically  done  in  meteorology.  The  formula  used  is 
identical  to  that  used  in  Smedstad  et  al.  (2003): 


A  C(f,x) 


£(/-/)(*-*} 

\fUf-T)2\/Ux-t)2 


where  /  is  the  model  mean  and  x  is  the  mean  of  the 
analysis.  In  this  case,  the  SSH  mean  of  the  model  and 


the  analysis  is  the  model  mean  over  1997  to  2002  with 
the  modifications  described  in  Section  2.3. 

Fig.  12  shows  AC  for  the  world  ocean  and  11 
subregions  averaged  over  22  30-day  forecasts  from  the 
1/16°  and  1/32°  models.  The  higher  resolution  model 
has  more  larger  amplitude  small-scale  features  which 
can  get  out  of  phase  on  shorter  time-scales.  Therefore, 
persistence  has  less  skill  in  the  1/32°  system  compared 
to  the  results  at  1/16°  resolution.  Thus,  the  best  way  to 
judge  performance  is  not  by  the  anomaly  correlation  for 
the  forecast,  but  by  the  spread  between  the  forecast  and 
persistence.  In  all  cases  the  spread  is  larger  for  the  1/32° 
model  on  any  day  where  the  difference  is  readily 
discemable  in  Fig.  12. 


Forecast  Length  (days) 


Fig.  12,  1/16°  (dashed  lines)  and  1/32°  (solid  lines)  SSH  forecast  verification  for  the  global  domain  and  1 1  subregions.  The  verification  statistic  is 
anomaly  correlation  between  the  model  forecast  (red)  or  persistence  (blue)  and  the  veil  tying  analysis.  The  plots  show  the  mean  statistics  over  22  30- 
day  forecasts  initialized  during  the  period  1  June  2001-31  May  2002.  See  Table  2  for  regional  boundaries. 
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Table  2 

Regional  boundaries  of  the  subregions  depicted  in  Figs.  12  and  14 


Subregion  name 

Latitude  range 

Longitude  range 

Gulf  Stream  region 

35°N-4$aN 

76°W-40eW 

Agulhas  retro  flection 

45°S-25&S 

0°-5G°E 

Equatorial  Atlantic 

20°S-20*N 

5S°W-i5°E 

Mediterranean  Sea 

30&N-47°N 

6°W-36°E 

Kuroshio 

25°N-50°N 

127°E- 180° 

NW  Pacific 

40°N-58ffN 

150°E-I80* 

Equatorial  Pacific 

20°S-20&N 

109°E-77°W 

Hawaii 

1&°N-235N 

159,5aW-t53,5°W 

Japan/East  Sea 

34QlSf-4g°N 

127*E-t42,5°E 

NW  Arabian  Sea  and 
Gulf  of  Oman 

15°N~26°N 

5 1*E-65*E 

Western  Arabian  Sea 
and  Gulf  of  Aden 

EQ-200N 

40*E-60*E 

To  further  compare  the  forecast  skill  of  the  1/1 6°  and 
1/32°  systems,  23  30-day  forecasts  of  SSH  from  each 
model  (initialized  during  the  period  l  June  01-31  May 
02)  were  compared  with  tide  gauge  time  series  at  78 
locations  (Fig.  13).  Tide  gauge  SSH  is  not  assimilated 
by  either  system.  A  1 3-day  moving  average  was  applied 
to  filter  time-scales  not  resolved  by  the  assimilated 
altimeter  data.  The  forecast  SSH  from  the  1/32°  system 
generally  agrees  better  with  the  tide  gauge  observations 
for  the  forecast  duration  at  both  open  ocean  (island)  and 
coastal  tide  gauge  locations. 

3.5.3.  Sensitivity  of  ocean  forecast  skill  to  the  forcing 
fields 

A  significant  question,  addressed  in  Smedstad  et  al. 
(2003),  is  the  sensitivity  of  ocean  model  forecast  skill  to 
the  use  of  climatological  atmospheric  forcing  fields 
beyond  the  end  of  the  forecast  available  from  an 
atmospheric  model.  This  question  is  revisited  here  in 
Fig.  14.  These  results  indicate  that  globally  and  in  many 
regions  of  the  world  ocean  the  30-day  SSH  forecasts  are 
not  very  sensitive  to  the  difference  between  the  clima¬ 
tological  and  analysis  quality  forcing  fields.  The  evo¬ 
lution  is  more  sensitive  to  the  initial  state  than  to  the 
atmospheric  forcing  anomalies  on  the  time-scale  of  the 
30-day  forecasts. 

Hurlburt  (1984)  discusses  classes  of  oceanic  response 
to  atmospheric  forcing  and  the  corresponding  implica¬ 
tions  for  the  time-scale  of  forecast  skill.  For  example,  the 
evolution  of  mesoscale  variability,  predominantly  due  to 
flow  instabilities,  typically  has  low  sensitivity  to 
atmospheric  forcing  on  30-day  time-scales.  Examples 
are  mesoscale  eddies  and  the  meandering  and  eddy 
shedding  of  ocean  currents  and  fronts.  These  are 
strongest  in  regions  such  as  the  Gulf  Stream,  Kuroshio 
and  Agulhas  retro  flection,  but  are  ubiquitous  in  the 


world  ocean  (e.g.,  see  Figs.  2,3  in  the  northern  Arabian 
Sea  and  Fig.  9  in  the  Hawaiian  Islands  region).  In 
equatorial  regions,  shelf  break  wave  guides,  and  in 
shallow  water  (the  latter  mostly  excluded  in  NLOM),  the 
model  response  to  the  wind  forcing  is  more  deterministic 
and  the  SSH  forecast  skill  improves  when  analysis 


Fig.  13.  A  comparison  of  forecast  and  observed  SSH  fluctuations  from 
the  i/16*  (blue  lines}  and  1/32°  (red  lines}  NLOM  systems,  A  13-day 
moving  average  was  applied  to  filter  time-scales  noi  resolved  by  the 
altimeter  data.  The  median  correlation  statistics  (computed  over  23 
forecasts  initialized  between  l  June  2001  and  31  May  2002)  are 
computed  over  (A)  tide  gauges  at  open  ocean  island  stations  (49  total), 
(B)  coastal  tide  gauges  (29  total)  and  (C)  all  tide  gauges  (78  total).  Tide 
gauge  SSH  data  are  not  assimilated  by  the  systems 
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Fig.  14.  $SH  anomaly  correlation  of  1/32°  NLOM  30-day  forecasts  with  the  verifying  analysis  for  the  global  domain  and  M  subregions.  The  red 
carve  is  for  the  NLOM  forecast  made  with  operational  forcing;  the  green  curves  signify  the  forecasts  which  used  analysis  quality  wind  and  thermal 
forcing  and  SSI  fields,  and  the  cyan  curve  verifies  forecasts  of  persistence  (Le.,  no  change  from  the  initial  state).  The  results  arc  means  for  26 
forecasts  initialized  during  the  time  window  l  June  2001-31  May  2002. 


quality  winds  are  used  for  the  duration  of  the  forecast. 
Even  in  this  ease,  the  propagation  of  ex  isting  equatorial 
and  coastal  ly- trapped  waves  can  be  forecast  with  some 
skill  after  the  reversion  to  climatological  forcing,  Jacobs 
et  al.  (1994)  even  demonstrated  —  7-year  skill  when 
climatological  forcing  was  used  while  predicting  the 
propagation  of  a  trans-Pacific  Rossby  wave  (along  30°N) 
that  was  generated  by  the  1982-83  El  Nino, 

This  is  a  more  definitive  investigation  of  the 
sensitivity  of  forecast  skill  to  atmospheric  forcing  than 
Smedstad  et  aL  (2003),  who  used  8  forecasts  early  in 
2001  and  fewer  subregions  to  address  this  question.  In 


this  comparison  a  total  of  26  forecasts  were  used  with 
initial  states  covering  a  complete  annual  cycle  from  June 
2001 -May  2002.  Also,  the  quantity  and  quality  of  the 
altimeter  data  are  significantly  improved* 

4,  Summary  and  conclusions 

A  1/32°  nearly  global  ocean  nowcast/forecast  system 
has  been  developed  by  the  Naval  Research  Laboratory 
at  Stennis  Space  Center.  This  system  uses  the  NRL 
Layered  Ocean  Model  (NLOM)  with  assimilation  of  sea 
surface  height  from  satellite  altimeters  and  sea  surface 
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temperature  from  multi-channel  satellite  infrared  radio¬ 
meters.  This  system  has  been  running  in  near  real-time 
since  1  Nov.  2003  and  daily  in  real-time  since  l  Mar. 
2005.  It  became  an  operational  system  at  the  Naval 
Oceanographic  Office  on  6  March  2006,  replacing  the 
existing  1/16°  system  which  ceased  operation  on  12 
March  2006.  Real-time  and  archived  results  are 
available  online  at  http://www.ocean.nrlsscmavy.mil/ 
globaLnlom. 

The  1/32°  system  has  improvements  over  the  earlier 
system  that  can  be  grouped  into  two  categories:  better 
resolution  and  representation  of  dynamical  processes 
and  design  modifications.  The  design  modifications  are 
the  result  of  accrued  knowledge  since  the  development 
of  the  earlier  1/16°  system.  Some  of  these  modifications 
could  have  been  implemented  in  the  1/16°  system,  but 
were  included  only  in  the  1/32°  system  because 
updating  the  1/16°  system  and  developing  the  1/32° 
system  concurrently  wasn't  feasible  given  manpower 
and  computational  constraints. 

The  improved  horizontal  resolution  of  the  1/32° 
system  has  important  dynamical  benefits  which  have 
implications  for  the  model's  ability  to  accurately  nowcast 
and  skillfully  forecast.  At  the  finer  resolution,  current 
pathways  and  their  transport  become  more  accurate,  the 
global  ocean  circulation  is  changed  (including  inter-basin 
exchange)  and  the  surface  SSH  variability  increases  and 
becomes  more  realistic.  These  improvements  result  in  the 
1/32°  system  being  a  better  dynamical  interpolator  of 
assimilated  satellite  altimeter  track  data,  using  a  1-day 
model  forecast  as  the  first  guess  for  the  assimilation.  The 
result  is  quantitatively  more  accurate  nowcasts,  as  is 
illustrated  by  the  model -data  comparisons  presented  in 
Section  3.  Based  on  comparisons  with  ocean  color 
imagery  in  the  northwestern  Arabian  Sea  and  the  Gulf 
of  Oman,  the  1/32°  system  even  demonstrated  the 
capability  to  map  small  eddies,  25-75  km  in  diameter, 
with  70%  reliability  and  a  median  eddy  center  location 
error  of  22.5  km,  a  surprising  and  unanticipated  result 
from  the  assimilation  of  altimeter  track  data.  For  all  the 
eddies  (50%  small  eddies),  the  reliability  was  80%  and  the 
median  eddy  center  location  error  was  29  km.  The  1/1 6° 
system  was  successful  for  eddies  >  50  km  in  diameter  and 
gave  a  median  eddy  center  location  error  of  32.5  km 
for  the  small  eddies  it  depicted.  The  smallest  eddies 
the  models  could  realistically  represent  were  8><6 
grid  intervals  in  diameter  or  1/4°  for  the  1/32°  model 
and  1/2°  for  the  1/16°  model  at  20°N,  near  the  middle 
latitude  of  the  region  used  for  this  comparison. 

The  1/32°  system  exhibits  improved  forecast  skill  in 
relation  to  the  1/16°  system.  This  is  due  to  (a)  a  more 
accurate  initial  condition  from  which  to  start  the  forecast 


and  (b)  better  resolution  and  representation  of  critical 
dynamical  processes  (such  as  upper  ocean  -  topographic 
coupling  via  mesoscale  flow  instabilities)  which  allow  the 
model  to  more  accurately  evolve  these  features  in  time 
while  running  in  forecast  mode  using  forecast  atmospher¬ 
ic  forcing  for  the  first  5  days,  then  gradually  reverting 
toward  climatology  for  the  remainder  of  the  30-day 
forecast  period.  At  1/32°  resolution  forecast  SSH 
generally  compares  better  with  unassimilated  tide  gauge 
observations  and  the  anomaly  correlation  of  the  forecast 
SSH  shows  a  larger  spread  when  compared  with 
persistence  than  the  1/16°  system. 

Applications  for  this  system  include  assimilation  and 
synthesis  of  global  satellite  surface  data;  ocean  predic¬ 
tion;  optimum  track  ship  routing;  search  and  rescue; 
anti-submarine  warfare  and  surveillance;  tactical  plan¬ 
ning;  high  resolution  boundary  conditions  that  are 
essential  for  even  higher  resolution  coastal  models;  sea 
surface  temperature  for  long  range  weather  prediction; 
inputs  to  ice  models,  biochemical  models  and  shipboard 
environmental  products;  environmental  simulation  and 
synthetic  environments;  observing  system  simulation 
and  assessment;  ocean  research  and  education;  pollution 
and  tracer  tracking;  fisheries  and  marine  resource 
management;  impact  on  ocean  structures  such  as  oil 
rigs  and  inputs  to  water  quality  assessment. 

The  capability  to  map  small  eddies,  25-75  km  in 
diameter,  demonstrated  by  assimilation  of  altimeter  data 
into  the  1/32°  model  has  strong  implications  for  nested 
finer  resolution  coastal  region  models,  including  (I)  the 
value  of  altimeter  data  for  assimilation  by  these  models 
and  (2)  the  value  that  high  resolution  boundary 
conditions  can  provide  them.  These  models  should 
give  enhanced  capability  in  mapping  small  eddies,  using 
altimeter  data  assimilation,  over  the  1  /32°  model  due  to 
finer  resolution  and  to  stronger  topographic  and 
atmospheric  forcing  constraints  than  found  in  the  open 
ocean.  However,  it  will  be  essential  to  separate  the 
temporally  unresolved  wind-driven  non-steric  compo¬ 
nent  of  sea  surface  height  from  the  altimeter  data, 
especially  in  shallow  water.  Doing  this  should  also 
benefit  altimeter  data  assimilation  into  global  ocean 
models.  The  separation  could  be  done  using  a  global 
high  resolution,  wind-driven  baro tropic  ocean  model 
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